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Abstract 



(N 

^ ' The aim of this paper is to present a governing equation for first order axial metric per- 

turbations of general, not necessarily static, spherically symmetric spacetimes. Under the non- 
^5 ' restrictive assumption of axisymmetric perturbations, the governing equation is shown to be a 

two-dimensional wave equation where the wave function serves as a twist potential for the axisym- 
metry generating Killing vector. This wave equation can be written in a form which is formally 
a very simple generalization of the Regge- Wheeler equation governing the axial perturbations of 
a Schwarzschild black hole, but in general the equation is accompanied by a source term related 
to matter perturbations. The case of a viscous fluid is studied in particular detail. 

U '. PACS: 4.30.Nk, 4.40.Dg, 4.40.Nr 



1 Introduction 

To consider perturbations of various physical fields around a spherically symmetric background is natural 
for several reasons. From the mathematical point of view, the symmetries of the background simplify the 
equations considerably and make them tractable for taking the analytical treatment further and for carrying 
out the numerical calculations needed without going far beyond standard methods. From the physical point 
of view, many systems (in this paper, the ones in mind are mainly astrophysical) owe their main structure 
to their monopole moment, and may consequently be regarded as perturbations of spherically symmetric 
systems to a good approximation. Moreover, as far as e.g. gravitational and electromagnetic fields are 
concerned, being non-radiative when restricted to be spherically symmetric, it is of course highly interesting 
to see what types of radiation are allowed for their more general neighbours. 

The discussion of this paper will be restricted to integer spin fields, satisfying linear differential equations 
defined on a spherically symmetric background. All such fields are known to split naturally into two parts 
that are not coupled by the equations of motion, namely an odd parity (or axial) part and an even parity 
(or polar) part. A standard way of separating out the angular dependence from the equations of motion 
for such fields is to make use of scalar, vector and tensor spherical harmonics (cf. ||l|), which all can be 
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defined in terms of the scalar ones, i.e. Yf^{d, (/>), and derivatives tliereof. As the final form of the equations 
are known not to depend on the separation constant m, there is no mathematical restriction of setting m 
to zero from the outset, corresponding to only considering axisymmetric fields. This was for instance the 
approach taken by Chandrasekhar in his extensive treatment of gravitational perturbations of Schwarzschild 
and Reisser- Nordstrom black holes Q , which summed up and added insight to the original works of Regge 
and Wheeler |^ and Zerilli Q|. Later, the same approach was taken by Chandrasekhar and Ferrari Q when 
generalizing to static spherically symmetric perfect fluid stellar models. The main focus of this paper is 
to use the axisymmetric approach to present a general treatment of the axial gravitational perturbations, 
valid for all spherically symmetric backgrounds and with no restriction on the types of matter perturbations. 
Only after this work was completed, this author became aware of the work of Gerlach and Sengupta ||^ who 
presented such a treatment for both polar and axial perturbations, using a 2 + 2 splitting scheme. However, 
new approaches and results are given in the present paper. Firstly, the method used to derive the final form 
of the perturbation equations, being a rather straightforward application of the general covariant form of 
the perturbed Einstein equations, differs from the one used in Here, the axial perturbation equations 
are first covariantly reduced to a Maxwell-like equation on the background geometry, valid not only for 
spherically symmetric backgrounds, but for any background with a non-null and twist-free Killing vector that 
is perturbed in such a way that the Killing vector obtains a twist to first order. The Killing vector is in this 
work identified with the axisymmetry generator of the neighbour to the spherically symmetric background, 
but it could e.g. just as well be identified with the time-like Killing vector of a stationary neighbour to a 
static background, meaning that the Maxwell-like equation can also be used in other contexts. Secondly, 
the further reduction of perturbation equations is shown to always lead to a governing wave equation of the 
Regge- Wheeler type {—di — V)'ip — S, with t and x being harmonically conjugate coordinates on the 2- 
space orthogonal to the SO (3) group orbits, V being a potential only depending on the background geometry 
and S being a source term directly related to the perturbation of the stress-energy tensor. It is found that 
the potential V can be written in a form which is a strikingly simple generalization of the form found by 
Chandrasekhar and Ferrari for the case of static perfect fluids: the energy density and pressure combination 
p — p simply has to be replaced by p — p± , where p± is the radial pressure of a radially moving observer 
who measures p as the energy density. It is also made clear that the wave function tp simply corresponds to 
a two-dimensional twist potential of the axisymmetry generator, given of course that the perturbations are 
axisymmetric. Thirdly, for static perfect fluid backgrounds it is known that the source term S has to vanish 
unless the perturbation is stationary, in which case it corresponds to adding a small rotation to the fluid 
ball. Since it thus would be interesting to see what kinds of source terms are allowed when considering more 
general fluids, we discuss in detail a fluid with viscosity and find that S is in general non- vanishing, which 
corresponds to an interaction between shear and gravitational oscillations. 

In addition to the treatment of axial gravitational perturbations, we also consider Maxwell and Klein- 
Gordon test fields on a general spherically symmetric background. Also in these cases the equations of motion 
are found to be reducible to Regge- Wheeler type wave equations of the form given above. In fact, it is found 
that the potential V for the three types of fields can be written in the common form V = Vg + Vi, where 
Vs and Vi depends quadratically on the spin s and the multipole moment /, respectively. The Vi part can 
notably be identified with the centrifugal potential for null geodesies. Such a splitting of the wave potential 
into a spin part and a centrifugal (or angular momentum) part has earlier been discovered for Scharzschild 
black holes 0, but to the author's knowledge it has not been shown to generalize to the general spherically 
symmetric case. 

2 Spacetimes admitting a nearly twist-free Killing vector 

If a spacetime {M,gah) admits a Killing vector which is close to being twist-free, i.e. hypersurface or- 
thogonal, it is natural to treat it is as a perturbation of a spacetime admitting a Killing vector which is 
exactly hypersurface orthogonal. We shall therefore consider a one-parameter family of spacetimes with an 
associated family of Killing vectors 77°. The metrics of these spacetimes will be written in the form 

gab = J-afc + FpaP-b (1) 
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where 



F = ri'^Tla IJia = F-^-na (2) 

and ±a'' is the projection operator that projects onto the tangent subspaces orthogonal to 77°. Before 
continuing with deriving perturbation equations, it wiU be useful to write down some formulas originating 
directly from the Killing vector equation W^aVb) = in the case of 77" being non-null (i.e. F 7^ 0). To begin 
with, we note that the two-form ^aVb can be decomposed as 

Va?7b = V[aVb] + Vab (3) 

where 

77%a = ?7%ab = Vf^ab) = 0. (4) 

Contracting eq. (||) by if' it follows that 

Va^-Va\ri\F\. (5) 
In terms of = F^^rja, eq. can be equivalently expressed as 

Vq^6 = ^i{aVb) + F^^Vab (6) 

which immediately implies that Vab is proportional to the exterior derivative of /j.^; 

Vab = \FVlab ^ab = 2V[a/l6]. (7) 

Hence, we conclude that VaJyb is completely determined by F, VaF and the exact two-form U,ab as 

Va77fc = ~V[a^b] In \F\ + ^F^ab- (8) 

It should be noted that F and /ia, and hence their exterior derivatives VaF and fiat, are Lie dragged by 77°. 
Indeed, 

Cr,F = ?7'^Va(?7S&) = 27]" r]''V^aVb) = = F-^ C,r]a = (9) 

where Crii]a = of course follows from £,,77" — and -C.g"'' = 0. Also, it clearly follows from eq. (Q) that 
the projection operator _La'' (with any index positioning) is also annihilated by Now, it is easily realized 
that the Killing vector 77" is hypersurface orthogonal (aligned with a gradient) precisely when the two-form 
flab vanishes. Indeed, the failure of hypersurface orthogonality can in four dimensions be described by the 
twist vector ^ (sometimes defined with opposite sign as in j^) 

which according to eq. (^) is dually related to flab as 

17- = ^^Fe'^^'^Sb^.d ^ab = F-^abcd^'v''- (11) 

Consequently, fJ" also satisfies C,^fl°- — 0. 

Let us now consider a one-parameter family of axisymmetric solutions to the Einstein equations, with 
the aim of linearizing these equations around a particular solution for which flab = 0. Referring to A as the 
parameter of the family, the notation 6 = d/dX will be used for derivation of tensor fields. The first order 
metric purturbation tensor ^ab — Sgab turns out to split naturally into two parts according to 

lab = ^7a6 + ^lab 

+ -fab = (-La'' -Lb'' + fiaf^bV'^V'^hcd = S±ab + {SF)fiaf^b (12) 
"7ab = 2-q''^H^a^b)'^lcd = 2F/i(a(5/ib) = 277(^(5^;,). 
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As always, in the absence of a natural identification of spacetime points along the family of spacetimes (i.e. 
a preferred one-parameter family of diffeomorphisms), S is only defined up to a gauge transformation 



S^5 + C(; (13) 

given in terms of an arbitrary vector field However, in the present case it is natural to only consider 
the families of diffeomorphisms whose associated tangent space mappings takes the Killing vector 77° into 
"itself" , corresponding to the condition 

Sr]" = (14) 
which clearly implies that we only allow for gauge transformations with vector fields restricted to satisfy 

Lev" = -A,C° = 0. (15) 
It then follows from 7/°/ia = 1, V°'-Lab — and g'^'^gcb — S'^b that 

v''S^ia = (16) 

V"S±ab = (17) 
gac gbd Sl""^ ^ -Slab - 277(a(5/ih) (18) 

which will simplify the calculations to come considerably. To derive what we shall refer to as the axial 
perturbation equations, we split the Einstein equations in the form Zai, '■= Rab — i^{Tab — \T'^cgab) = into 
its "components" orthogonal to 77° according to 

l_ac^bd2^^ = (19) 

ri'^ri^Zab - (20) 
±'^''7,^Zb, = 0. (21) 

As we are now about to apply 5 to these equations, we shall need the general form of SRab, which reads 

5Rab = -iVaVfc7^c - iQTab + (alb)c (22) 

where □ ~ VVq. Now, it does not take a too frustrating amount of algebra to realize that when evaluating 
eqs. ( |r^ - (^) at ^lab = 0, using the above derived formulae, the ^^ab part of the metric perturbation tensor 
7afc satisfies eqs. (^9|) and (|o|) identically, while on the other hand ^7ah satifies eq. (^ ) identically. Hence 
it follows that the perturbation equations separate ^^ab and ^7a6, implying that those parts of the metric 
perturbation can be treated separately. We shall here only focus on the ^^ab part, later to be identified 
with the axial (or odd parity) part when we specialize to the case of a sperically symmetric background. We 
hence insert ^ab — ~Jab into eq. (^), evaluate at flab = to find the Maxwell-like equation 

Vfa(FQ'^'') = W (23) 

where 

Qab = Snab (24) 

ja ^ 2ST'' T" = l-^^r^^Tbc. (25) 

It should be noted that since flab and vanish on the background - the former field by assumption, the 
latter by the twist-free Killing vector ry" being an eigenvector of the Ricci tensor and hence also of the 
stress-energy tensor (cf. 0) - it directly follows from eq. ( p^ ) that the perturbation fields Qab and J" are 
invariant under the gauge transformation (p^). This can be recognised as a lemma by Stewart and Walker 
pO{ : the linear perturbation of a quantity is (identification) gauge invariant if that quantity vanishes when 
unperturbed. Moreover, using (5r/° = it can easily be shown that Qab and J° share the properties of 
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their respective nonlinearized fields of being orthogonal to, as well as Lie dragged by, 77°. The perturbation 
equation ( ]23| ) should be accompanied by these conditions, as well as the condition of Qab being closed, which 
follows from the closedness of Q,ab and the fact that 5 commutes with exterior differentiation. Moreover, as 
an integrability condition of eq. (p3[), it follows that the linearized matter current J" is conserved, 

VaJ" = 0. (26) 



In the next section we will in fact not be working with the Maxwell-like equation (|2^) directly, but rather 
with its dual equation 

2V[aQ6] = -nJab (27) 

where 

= Sn'' Jab = eabcdJ^v''- (28) 

The linearized twist vector as well as the dual matter field Jab can easily be shown to also have the 
properties of being orthogonal to and Lie dragged by 77°. It should be noted that the duality relations 
( |rT| ) for Vlab and fi" holds also for their linearized versions Qab and which is due to the fact that the 
perturbation equations are evaluated on a background where ^ab and 17" vanishes. Moreover, the above 
relation between Jab and J° can be inverted into 

ja ^ _lp-l^abcd^^j^^^ (29) 

Finally, in this dual formulation, the additional equations '^[aQbc] = and VaJ"' — are replaced by 

VaiF-^Q") - (30) 

^[aJbc] = 0. (31) 



3 Axial gravitational perturbations 

We shall now restrict the reduced perturbation equations derived in the previous section to the case of a 
spherically symmetric background perturbed into an axisymmetric neighbour. As mentioned in the intro- 
duction, the Killing vector 77° will thus be chosen as the generator of the axisymmetry, reducing simply to an 
SO (3) generator for the background. The metric of the background geometry will be written in the standard 
form 

ds^ = da^ + (de^ + si-a^e dc/)^) (32) 

where da^ is the induced Lorentzian two-metric on the space orthogonal to the SO (3) group orbits and r is 
the Schwarzschild radius. For simplicity, we will arrange the coordinates so that the Killing vector 77° will 
coincide with {d/d(f)°'. Note that the scalar F = ri°"qa on this background takes the form 

= (r sin 61)2. (33) 

To reduce the perturbation equations to two-dimensional ones by separating out the dependence of the 
remaining angular variable 6, we start off by noting that no restricion is implied by assuming that the two- 
form Jab, which by eq. (|3^) is closed, can be given in terms of a potential Ya which is orthogonal to the 
SO (3) orbits (i.e. whose 9 component is vanishing). Thus we set 

Jafc = 2V[,n]. (34) 

It then follows from eq. ( p7| ) that the linearized twist vector can be written in terms of Ya and a scalar 
<& as 

Qa = - KYa. (35) 
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Turning to eq. ( pO[ ) which corresponds to the closedness condition for Qab, the 9 dependence can be separated 
out by setting 

$ = C{e)rib Ya = C{e)Xa (36) 

where the scalar -0 and the one- form Xa have no angular dependence, i.e. they are Lie dragged by all SO (3) 
generators. Inserting eq. (|3^) into eq. ( ^o|) leads to a separation into an ordinary differential equation for 
C{0) as well a two-dimensional wave equation for -if) coupled to Xa- The equation for C{9) will be discussed 
further in the appendix. Here we merely note that it is solved by setting 

c{e) ^ G;^i\cose) (37) 

with Gi_^2^{y) being the ultraspherical (or Gegenbauer) polynomial. The wave equation can be covariantly 
written in form 

{VaV' - U)ij = kS (38) 

where Da is the connection of the Lorentzian two- metric da^. The potential U can be put in the formally 
simple form 

1 6m, 1(1 + 1) , , 

^=¥^- — ^^^ (39) 

with T invariantly defined (given the spherical symmetry) by ^kt being the eigenvalue function of the Ricci 
tensor with respect to the eigenvector rf^ (or any other S0(3) generator), i.e. 

R'^i.rf' ^Iktt^". (40) 

The reason for defining t with the factor inserted as above is that the background Einstein equations 
imply that r is simply minus the trace of the 2x2 block of the stress-energy tensor orthogonal to the 
SO (3) orbits, an interpretation which has a clearer physical significance. The function m is the spherically 
symmetric mass function JTTI , invariantly defined in terms of the Schwarzschild radius r according to 

2rT? 

VarW^l . (41) 

r 

The source term S on the left hand side of eq. (^) is given by 

S^rVa{r~^X'') (42) 

and is thus directly related to the perturbation of the stress-energy tensor. We will have more to say about 
the wave equation ( |38| ) , but this will be postponed to section ^ 

Having found the general solution to the conserved current condition Vq J" = by introducing a gauge 
potential Yq for the dual field Jab according to eq. (^4|), it is of interest to transform this solution back to 
J° which is more directly related to the matter perturbations. Using eqs. (|2^), ( ^4| ) and (^6|), it follows that 
J" can be expressed as 

r = (r^ sine*)-! [c'{e)j'' - c{e)j [d/dey] (43) 

where 

J» = e-«'X6 ^ Xa = eabJ'' J = e°'''VaXb (44) 

with tab being the Levi-Civita volume form of the induced two- metric. It directly follows that Va J'^ = can 
be solved by putting J"^ in the form given by eq. (^) and requiring that J'^ and J satisfy 

VaT = J (45) 
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as implied be eqs. (Q. It should be stressed that since the form J° may be restricted by the types of matter 
studied, it may not be practical to solve the conserved current condition by expressing JZ" and J in terms 
some Xa according to eqs. (|4^). Rather, the two-dimensional equation (^) in general has to be considered 
as nontrivial equations of motion for the perturbed matter fields. Hence it is appropriate to note that the 
source term of the wave equation ( ^8[ ) can be expressed in terms of J"" rather than as 

^-re-^'P^r-'jb). (46) 

Let us now discuss Israel's junction conditions of the derived perturbation equations. Returning 
first to the general, not necessarily spherically symmetric case, we assume that we are given a non-null 
hypersurface E with unit normal n", n^n^ = e = ±1, which may be given an arbitrary smooth continuation 
off of E. For practical purposes we need only consider the case when rf is tangent to the hypersurface, i.e. 
rfua = 0. The fields to be matched across S are the first and second fundamental forms defined as 

jab = gab - e riaUb (47) 
kab = ja'jb''^cnd. (48) 

By using a gauge transformation generated by a vector field C" satisfying 

V^Ca = CnC = (49) 

we may set Sria to zero. Such a gauge transformation only affects the polar part of the perturbations, since 
it leaves the axial field Sfia — S{F~^r]a) invariant. Hence we are able to make the gauge choice Sua = 
without making any gauge fixing in the axial sector that we are considering. For simplicity, we therefore set 
Sua = to find that 

Sjab = Sgab = lab- (50) 

The relevant projection for the axial perturbation junction condition is 

^a^^'^ib, = FSlla (51) 

from which we conclude that all components oiSfia should be matched. Note, however, that Sjia is not a gauge 
invariant field, since Ha is not zero on the unperturbed spacetime. Namely, under a gauge transformation 
generated by a vector field C° such that 

C = /?7" V'^^af = (52) 

one finds that (5/ia transforms as 

S^Ia + C^Ha = S^Ia + /A,Ma + (?7^M6)Vq/ = S^a + ^af- (53) 

As far as the second fundamental form is concerned, the relevant quantity to consider and match across E is 

5{^a''il'hc) = -\FQabn^ = -\F-^eabcdn''QW'- (54) 

Let us now specify these matching conditions to the spherically symmetric case. Excluding the nonradiative 
case I = 1 which will not be treated here, it is implied by eqs. (Bq) - (521) that J/i^ and Qab take the forms 



^A^a = - ,, , ^..f'^^i , . . eJ'[V,{ri^)-KX,] + Vaf (55) 
(( + 2)[l — l)r^ siYi 



Qab = eafc + 4^ (V[,0)e,f [2?,(r^) - kX,] (56) 

sm sm y 

where the term Va/ in the expression for 5[ia corresponds to the freedom to make the gauge transformation 
(m). Now, since (5/ia is to be continuous and since Vq/ — Vaf + df/dOVa^, it directly follows that df/d9 
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has to continuous. Furthermore, the possibihty of adding to / a ^-independent term with discontinuous 
gradient is excluded by the fact that C'{9)/ s\v?9 is not constant for I ^ 1. Consequently, the matching of 
Sjia leads to the condition that eo'^\Vb{ril}) — nXb] and hence 'Da{r%li) — KXa be matched. For all practical 
purposes we can assume that n° is orthogonal to the S0(3) orbits, implying that n°Va0 = and hence 

Qabu' = Cabu' + ^b'[Dc{rij) - nX^^. (57) 

r-^ sm 6' sm 

It directly follows that the matching of Qatn'' leads to the further condition that the two-scalar ip itself 
should be matched. 



4 Klein- Gordon and Maxwell fields 

In this section, we shall consider Klein-Gordon and Maxwell test fields on a general spherically symmetric 
background and, as was done for the axial gravitational perturbations treated in the previous section, reduce 
the equations of motion to two-dimensional wave equations. 

Starting with the Klein-Gordon case, the equations of motion is taken to be 

□$ = eoj (58) 

where an unspecified source scalar j multiplied by a coupling constant Eq has been added to the free massless 
Klein-Gordon equation. Restricting to axisymmetric fields, $ as well as j will be assumed to be Lie dragged 
by r]", one of the S0(3) generators. Referring to the form (32) for the general spherically symmetric metric, 
the 9 dependence can straightforwardly be separated out according to 

$ = C(6l)r-V (59) 
J = C{9)r-'S (60) 

with tp and S being S0(3) invariant scalars. Taking C{9) to be the Legendre polynomial P;(cos6'), the 
Klein-Gordon equation reduces to the covariant two-dimensional equation 

{VaV" - U)^ = eoS (61) 

which is formally identical to the eq. (p8|), but with the first two terms in the potential U having different 
coefficients; 

1 2m 1(1 + 1) 

As for Maxwells equations 

^[aFbc] = (64) 
Vaf - (65) 

with the electromagnetic field strength Fab and the current j° assumed to be Lie dragged by 77", the hyper- 
surface orthogonality of 77" makes it natural to make the splitting 

Fab = Pab + 2M[QS'f,] (66) 
ja = Ja + JfJ-a (67) 

where all introduced fields are assumed to be orthogonal to as well as annihilated by Cr/- Indeed, using 
S7[afJ-b] — leads to a decoupling into the two sets of equations 

FVa{F-'S'')^eiJ (68) 
V[aSb] = (69) 
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and 

VbP"'' ^ eiJ" (70) 

V[aPbc] = (71) 

VaJ" = (72) 

corresponding to the axial and polar sector, respectively. 

Starting with the simpler axial case, we solve eq. ( |69| ) by setting 

Sa = Va$. (73) 

The angular dependence can now be immediately be separated out, this time according to 

$ = c{e)ij (74) 

J^C{9)S (75) 

with C{9) taken to be the ultraspherical polynomial G^^l (cos 9) and ip and S having no angular dependence. 
The resulting two-dimensional equation again takes the form 



{VaV'' - U)^ = ei^ (76) 
but with the potential U in this case only consisting of the centrifugal term 

U= ^ \ (77) 

Turning to the polar case, which clearly can be treated similarly to the axial gravitational perturbations, we 
introduce the fields P° and Jab dual to Pab and J° according to 

1 ^abcd^ 



pa ^ _i^abcct^^p^^^ p^^ ^ F-^abcdP'v" (78) 

Jab = eabcdJV J" = -\F~^e''^^SbJcd- (79) 

This turns eqs. ^ - (|7|) into 

2V[„P6] = -ei Jab (80) 

VaiP-^P") = (81) 

^aJbc] = 0. (82) 

We start off by solving the closedness condition for Jab according to 

Jab = 2V[,n] (83) 

with the condition that the gauge potential Ya be orthogonal to the SO (3) orbits. According to eq. (|80|), the 
field Pa can be written as 

Pa = Va$ - eiYa (84) 

for some scalar potential <&. To reduce the remaing equation (|8l| ) into a two-dimensional one, we set 

$ = C{9)i^ Ya = C{9)Xa (85) 

with C{9) = (cos 6*) and with if) and Xa being S0(3) invariant fields. Again, we obtain a two- 

dimensional wave equation of the form (p8|). As for the axial Maxwell case, the potential U is given by eq. 
([fz]), while the source function S has the form 

5 = VaX"- = e-^'^VaJb. (86) 
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The field J a here refers to the form of J°, 

r = (r^ sin6')-i [C"(6') J° - C[B)J (8/89)"] (87) 

where, as implied by the duality between J° and Jab, 

J- = e-^Xb J = e'^'Pa^b. (88) 

In analogy with the axial gravitational perturbations, the conserved current condition Va J" is thus found 
to reduce to the two-dimensional equation 

VaJ" = J. (89) 

Having found that the linear equations of motion for Klein-Gordon fields. Maxwell fields (axial as well 
as polar) and axial gravitational perturbations are all governed by a two-dimensional wave equation of the 
form 

(P"Pa - U)^ = eS (90) 

it is interesting to note that the potential U can be naturally split into a spin dependent part Us and an 
angular momentum dependent part Ui as 

U = Us + Ui (91) 



Us^{s- 1) 



1 2m 

- (s + 1)^ 



(92) 



^ ^ (93) 

with s taking the values 0, 1 and 2 for scalar, electromagnetic and linearized gravity fields, respectively. 
Note also that the value of the coupling constant e on the left hand side of eq. ( p0| ) depends on the spin s 
as well; e = Ss with £2 = k. In geometrical units 2ei = £2 = 87r. Now, to connect more closely to previous 
works on various static backgrounds, let us choose as coordinates two harmonically conjugate functions |^ 
t and X to make the two-metric da^ take the manifestly conformally flat form 

da^ = e^''{-dt^ + dx^). (94) 

This turns the wave equation ( po|) into the form 

&^ 

dt'^ dx 

where the potential V and the source H are simply given by rescaling the invariant potential U and source 
S with the conformal factor e^", i.e. 

V = e^^U H = e^^S (96) 

implying of course that V also naturally splits according toV — Vg + Vi. The coordinates t and x are clearly 
not uniquely defined, since the freedom of choosing them is in correspondance with the infinite dimensional 
conformal group in two dimensions, but in the static case there is a natural particular choice. Starting out 
from Schwarzschild coordinates 

da^ = ^e^'^df + (1 - ^) ^dr^ (97) 

with V and m being functions of r only, the harmonically conjugate variables can preferably be taken to be 
the static time t and the Regge- Wheeler (or tortoise) radial variable 



e-''{l~^2rn/r)-^/'^ dr (98) 
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which gives dcr^ = e'^'^{—dt^ + dri). In this case, we may turn eq. ( p5| ) into a one-dimensional Schrodinger 
equation (in general with a source term) by separating out the time dependence in the standard manner 
according to 

iP = e^"*Z(r,) H = e"^'Y{r,) (99) 

which results in 

Clearly, besides the source Y, this is a standard one-dimensional Schrodinger equation with the squared 
frequency corresponding to the energy. 

The splitting of the potential V into a spin dependent and an angular momentum dependent part has 
previously been noted in the case of a Schwarzschild background (cf. [Q), for which the matter term r of 
course vanishes and the mass function m is a constant, but it does not appear to be have been known to 
generalize to the case of a general spherically symmetric background in the simple manner found here. Indeed, 
whether or not the background is static, the only formal modification of the potential V compared to the 
axial gravitational wave potential found by Chandrasekhar and Ferrari for static perfect fiuid backgrounds 
1^ , is the simplest possible: the matter combination p—p has to be replaced by t, minus the trace of the 2x2 
block of the stress-energy tensor that is S0(3) orbit orthogonal. Clearly, r can be expressed as p — p±, with 
p and p± being the energy density and radial pressure measured by an arbitrary radially moving observer. 
It should be noted, however, that whereas the combination p ~ p± is invariant under a change of such an 
observer, the same does in general not hold true for the quantities p and p± themselves. 



5 An application to fluids with viscosity 

As a concrete and physically relevant example, we shall here consider the general Regge- Wheeler type wave 
equation for axial gravitational perturbations of spherically symmetric backgrounds, in the case when the 
matter can be described as a fluid with viscosity in the standard relativistic manner. The stress energy 
tensor for a viscous fluid is ||l^ 

Tab = pUaUb + (p- C^)hab ^ "iWab (101) 

where the various flelds involved are the fluid four- velocity u", the energy density p, the pressure p, the 
coefficient of bulk viscosity the coefficient of dynamic (or shear) viscosity 77, the projection operator 
hab = UaUb + gab, the expansion Q = 'Vau"' and the shear tensor aab — {'^cU(^a)hb)'^ — \Qhab- For the 
background solution, the four- velocity u°' will be assumed to be orthogonal to ry°, meaning that the matter 
is not rotating around the symmetry axis. In particular, in the spherically symmetric case the fluid will 
be assumed to be radially moving, thus automatically satisfying rj°'Ua — 0. We now note that the vector 
T°- = ^''''ri'^Tbc takes the form 

T"^ = (p + p - Ce) v'^u, ±'''ub - 2?7±'^S'<^6c. (102) 

Applying the perturbation operator S to this vector, we arrive at the following form of the linearized matter 
current J": 



2 



where 



J'' = ST'' ^ {p + p- (Q - 2ria)P u'' + r] [FQ^^'u,, - Fh''^Vb{F-^ (3) - /3u"] (103) 

13 = v'Sua (104) 
a = F-^ri'^ifaab- (105) 
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The popping up of the two- form Qab in the expression for J° originates in the perturbation of the shear aab, 
being defined in terms of the covariant derivative of the four- velocity Ua- More precisely, we have used 

Si^aUb) = S/aSub + Uc ST^b S{±\±\rdevl = -jFQ'^b (106) 

where V^ab refers to the Christoffel symbols. As we now specialise to the case of a spherically symmetric 
background, the decomposition of J° in accordance with eq. (^) is achieved by decomposing the scalar 
T]°'Sua according to 

r]''6ua = sm~^eC'{9)f3 (107) 

where C{0), as in section ^, refers to the Gegenbauer polynomial Gi_^2^ {cos 0), while /3 is a scalar without 
angular dependence. The resulting expressions for the two-dimensional fields and are 

J7" = 2r^{p + p- Ce - 2rj ct) u° - 2r/ r [V' + r^n^Vbir^^ /S) + r ii p] n"" (108) 



J 



^2r]{n"'Va{riP) + [(/ + 2){l - I) + 2Kr^{p + p - QQ - 27^ a)] /3} 



(109) 



where ria = u^^ba, implying that n° is the (up to sign) unique unit space- like vector which is orthogonal to 
both u° as well as the S0(3) orbits. We have also introduced the notation ii for the norm of the acceleration 
, and used that u° — iin'^ . For simplicity, we shall now further restrict the background by requiring that 
it be static. Using Regge- Wheeler variables, the two-metric dcr^ and the unit vectors u° and n"' take on the 
forms 



^ e-''{d/dr.,Y 



(110) 
(111) 
(112) 



with the gravitational potential v being a funcion of only. The time dependence of the perturbation fields 
is straightforwardly separated out according to 



(113) 



The assumption of staticity implies that both the expansion Q as well as the shear Gab vanishes on the 
unperturbed spacetime. This means that both 8 and a will be set to zero in eqs. (108) and (109), which 
inserted into the matter current integrability condition (|45|) leads to the ordinary differential equation 



drl 



dr^ 



W B 



r 



where, after using the unperturbed Einstein equations, 

1) 



W 



K, , 2m 1(1 

o(^ + 5p) + -r + — 



2W1 



2m ii 
r r 



— e'ip + p). 



(114) 



(115) 



In obtaining eq. (114), wc have multiplied eq. ( p5| ) by the coefficient of shear viscosity tj. The case 77 = is 
hence recovered from eq. (114) by first multiplying by 77 and thereafter setting 77 = 0. This results in the 
simple condition 



LuB = 



(116) 



meaning that for vanishing coefficient of shear viscosity, the perturbation is either stationary (w = 0) or can 
be interpreted as a pure gravitational wave not interacting with the matter (B = 0). This is in agreement 
with Chandrasekhar and Ferrari's work on axial perturbations of isentropic perfect fluids B, but note that 
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we here have not assumed that the coefficient of bulk viscosity vanishes, so the fluid need not be an isentropic 
one for the same result to hold. 

Returning to the case ?? 7^ 0, the Schrodinger equation ( |100D (with s — 2), the source function Y can be 
written as 



Y 



-2r— [{p + p)e''B] + 2iujri 



dr. 



{r~^e''B) + e''Z 



(117) 



where it has been used that ii — e^^dv/dr^ to put Y into a compact form. Since Y contains the wave 
function Z , it is clearly appropriate to rewrite eq. ( |100D by moving the term in question to the left hand 
side, with the result 



-— -5- + V — uj{uj — 2iKrie'^) 
dri 



-kY 



(118) 



where 



V 



i^(p_p)__ + A_^ 



Y = -2r4- [{p + pVB] + 2iL0 7^r^e-'4-{'r~''^e''B). 
dr.. dr^ 



(119) 
(120) 



To summarize, the axial perturbations of a static sperically symmetric viscous fluid are governed by the 
system of two coupled ordinary differential equations ( |ll^ ) and (114). To solve the equations numerically, 
one starts at the center, using as initial conditions the following expansions in the Schwarzschild radius of 
the solution with regular center: 



Z = r'+'{ZQ + Z^r' 



B = r^+\Bo + Bir' + ...) 



(121) 
(122) 



where 



(2Z + 3)Zi 



;{2l - l)po -PO 



- 2'^o(wo - 2iKri) \ Zq 



k[{1 + l){po +po) - {I - l)iuJor]] Bo 



(123) 



{21 + '6)Bi = 



Wo 



1 + + 1) ) Po - [21 - l)po 



K 
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+ 7; — {pQ+Po)\ Bo 
I r\ ' 



(124) 



(125) 



To arrive at these expansion coefficients we have used 



P = Po + pir + 

„2 



m = ^Kpo + raxr"^ + . . . 
V = vo^ j^k(po + 3po)»'^ + 



(126) 



which follows from the background Einstein equations. The integration of the perturbation equations should 
be continued up to the stellar surface at r = i?, with R denoting the Schwarzschild radius of the star. While 
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the function B has no vacuum analogue, the integration of the function Z should be further continued into 
the vacuum region r > R where eq. (|ll§|) is replaced by the original Regge- Wheeler equation 



1 - 



2M 



6M 



l{l + l) 



Z = 



(127) 



where M = m\r=R. The general condition that -0 and Vair^p) — nXa be matched across r = R boils down to 

Z\r^R+ (128) 



A 

_d_ 

dr^ 



(rZ) + 2Ke''r^{p + p)B 



dr. 



(rZ) 



r^R+ 



^ 0. 



(129) 
(130) 



^R- 



Furthcrmore, if the condition of purely outgoing gravitational waves at infinity is imposed, corresponding to 
Z oc e"*'^'"* as r — > oo, one is lead to a discrete gravitational wave spectrum for the system, much in the same 
way as for perfect fluids stars. The difference is that the gravitational waves (closely related to the function 
Z) in this case interacts with fluid shear waves (closely related to the function B) in the stellar interior. 



6 Concluding remarks 

So far, most of the works on gravitational perturbations of spherically symmetric spacetimes have dealt with 
black holes or static perfect fluid stellar models, with focus on the determination of the quasinormal modes 
or means to excite them through various astrophysical processes (cf. j^] for a review). Gravitational radiation 
emitted by collapsing systems has also been studied to some extent, notable examples being the work by 
Cunningham, Price and Moncrief on Oppenheimer-Snyder colUapse |14, |l^ and by Seidel and Moore on 



more general collapsing perfect fluids ]16| | . Since Choptuik's discovery of critical phenomena in gravitational 
collapsing systems |l^ and the subsequent development of the subject into a brand new branch of general 
relativity, an increasing number of different types of matter has been studied in the context of, in particular, 
spherically symmetric gravitational collapse. Since analysing perturbations of the critical solutions is a 
crucial part of the game, it is clearly desirable to have a general framework for perturbations of spacetimes 
with spherical symmetry, but without any other restriction on neither the properties of the background nor 
on the perturbations themselves. It is hoped that this work, being complementary to Gerlach and Sengupta's 
1^, provides such a framework. 

In this paper, we have put the axial gravitational perturbations on an equal footing with Klein-Gordon 
and Maxwell test fields through the unified two-dimensional wave equation presented in the previous section. 
However, the more complicated polar gravitational perturbations have not been discussed at all. Although 
some attempts have been made to reduce the polar equations into a neat form, we have so far not been able 
to take the equations any further than Gerlach and Sengupta. However, judging from the vacuum case - a 
Schwarzchild background perturbed into a more general vacuum neighbour - one may suspect that a certain 
illuminating reformulation of the system of two-dimensional polar equations may be possible in general. 
Indeed, as discovered by Chandrasekhar Q, the vacuum axial and polar equations are related by what was 
later to be recognized as a supersymmetry duality [T^ , more precisely meaning that the Regge- Wheeler and 
Zerilli potentials for the axial respectively polar modes are supersymmetric partners with the implication 
that it is only necessary to investigate the properties of the former, simpler one. This indicates that it could 
be possible that the polar equations for the general spherically symmetric case can be reduced to an equation 
formally identical to the generalized Regge- Wheeler equation found here for the axial case, combined with a 
set of equations directly related to the matter perturbations. Whether or not this is merely wishful thinking 
remains an issue for further investigation. 

We have looked into the case of viscous fiuid stellar models in some particular detail and derived a 
system of equations for the axial perturbations that can be directly applied for numerical studies. It would 
be interesting to explicitly calculate the w-modes for some realistic models to see how the spectrum depends 



14 



on the shear viscosity rj which enters the equations. It would also be of interest to do a similar analysis for 
other types of matter sources that are more general than perfect fluids, such as elastic materials as described 
in a general relativistic setting by Carter and Quintana [|l9j . 
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Appendix: differential equation for the angular function C{9) 



For all integer spin fields studied in this paper, the separation of the dependence on the angular variable 9 
leads to the ordinary differential equation 

^ sini-2« e-^+k sini-2« 9) C{9) = (131) 
d9 d9 J 

where s ~ 0, 1, 2 is the spin of the field and fc is a separation constant. Setting k = [l + s){l — s + 1), the 
general solution can be written down as 

C{e) = (1 - y2)V2 [CiP« (cos 9) + CsQf (cos 9)] (132) 

where Pi{y) and QHy) are the associated Legendre functions of, respectively, first and second kind. To 
obtain solutions that are regular at y = ±1, we must set C2 to zero and I to an integer (which can be taken 
to be nonnegative since (^ + s)(^ — s + 1) is invariant under / ^ — Z — 1). The solution then becomes the 
ultraspherical (or Gegenbauer) polynomial G\l^~^{y), which for the radiative values of I, i.e. / > s, is related 
to the Legendre polynomial Pi{y) as 

G]il-\y) = N^{1 - yy-^Piiy) (133) 
where Nf is a normalization factor given by 

N? = l 

(2? - 1) (134) 

For the relevant values s ~ 0, 1, 2, it may also be noted that the Legendre polynomial can be expressed in 
terms of the ultraspherical polynomial according to 

Gy\y) = Pi{y) 

^yGJ-^riv) = -Piiv) (135) 



(/ + 2){l - l)Gl^i\y) + 2y-GJ^i\y) = -(1 - y')Pi{y). 
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